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Abstract 

The discovery and characterization of informative intraspecific genetic marl<ers is fundamental for evolutionary and 
conservation genetics studies. Here, we used nuclear ribosomal ITS sequences to access intraspecific genetic di- 
versity in 23 species of the genus Passiflora L. Some degree of variation was detected in 21 of these. The Passiflora 
and Decaloba (DC.) Rchb. subgenera showed significant differences in the sizes of the two ITS regions and in GC 
content, which can be related to reproductive characteristics of species in these subgenera. Furthermore, clear geo- 
graphical patterns in the spatial distribution of sequence types were identified in six species. The results indicate that 
ITS may be a useful tool for the evaluation of intraspecific genetic variation in Passiflora. 
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Introduction 

Understanding the distribution of alleles throughout 
the geographic range of a species is fundamental for molec- 
ular evolutionists, allowing inferences about the influence 
of historical processes in the spatial distribution of particu- 
lar lineages (Emerson et al, 2001). Phylogeography is a 
young discipline that represents a bridge linking population 
genetics, molecular phylogenetics and biogeography, 
amongst other fields (A vise et al, 1987; Avise, 2000). 
Phylogenetic methods can be used to construct haplotype 
trees that would indicate the historical relationships of gene 
lineages in a population or species. By comparing such 
trees to the geographical structure of the data one can infer 
historical patterns of population subdivision, as well as 
come to an understanding of current distribution of the ma- 
jor genetic clusters in a studied species, which may be use- 
ful for conservation purposes (Frankham et al, 2002). 
Furthermore, well-resolved phytogenies at the species- 
level or below are necessary to identify the taxa and popula- 
tions which have undergone recent speciation and those 
which are currently diverging (Bradshaw et al, 1995; 
Whittall et al, 2006). 
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For over two decades, nuclear ribosomal ITS (inter- 
nal transcribed spacers of the large subunit of ribosomal 
DNA) has been the most popular molecular marker in the 
nuclear genome for evolutionary studies in various plant 
groups {e.g. Baldwin et al, 1995; Alvarez and Wendel, 
2003; Hughes et al, 2006; Nieto-Feliner and Rossello, 
2007). Even though the use of ITS as a inolecular marker 
might be problematic because of paralogy and other com- 
plex evolutionary patterns caused by its highly repetitive 
nature (Alvarez and Wendel, 2003), the general evolution- 
ary mechanism for this region proposes that most intra- 
genomic copies would share identical sequences as a result 
of concerted evolution, a homogenizing process that de- 
pends on unequal crossing over during meiosis and biased 
gene conversion (Dover, 1994; Liao, 1999). 

Another advantage in using ITS as molecular marker 
in plants is that it provides an alternative to cpDNA mark- 
ers, which may inform only part of the history of a species, 
since cpDNA is usually inherited from only one of the par- 
ents (Birky Jr., 2001 ; Xu, 2005), and thus proposed migra- 
tion patterns based entirely on cpDNA may be inaccurate 
when applied to the population as a whole. The ITS region 
has been successfully used to infer phylogeographic pat- 
terns in a wide range of species (Jeandroz et al. 1997; 
Manos et al, 1999; Rodriguez-Lanetty and Hoegh-Guld- 
berg, 2002; Duran et al, 2004; Lorenz-Lemke et al, 2005; 
Koehler-Santos et al, 2006; Nettel and Dodd, 2007; 
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Yamaji et ah, 2007). Furthennore, since the ITS region is 
flanked by well-conserved rRNA genes, universal primers 
can be used for widely different plant groups, thereby 
avoiding the need for developing specific primer sets, as is 
the case for SSR markers. 

Passiflora L. is the largest genus in the Passiflo- 
raceae, consisting of over 520 species, recently split into 
four subgenera (Feuillet and MacDougal, 2003). The mo- 
lecular phylogeny of the genus has been investigated by 
various authors, and in general is in agreement with the 
morphological proposition of subgenera classification 
(Muschnere/fl/., 2003; Yockteng andNadot, 2004; Hansen 
et al, 2006). Particularly, Muschner et al. (2003) and AK 
Hansen (PhD. Thesis, University of Texas, 2004) have used 
ITS to estimate phylogenies for the whole genus. Krosnick 
and Freudenstein (2005) have used ITS with success to ana- 
lyze Supersection Disemma (subgenus Decaloba (DC.) 
Rchb. Lorenz-Lemke et al. (2005) and Koehler-Santos et 



al. (2006) have also assessed intraspecific variability for 
two Passiflora species based on ITS sequences. 

In this study, we used ITS sequences to evaluate gen- 
eral patterns of intra and interspecific variation in several 
species representing all four Passiflora subgenera, to test 
the usefulness of this genetic marker for further phylo- 
geographic studies in this genus. We found general differ- 
ences in ITS structure between the two most speciose 
subgenera {Passiflora and Decaloba), besides identifying 
geographic associations among lineages in at least six spe- 
cies. 

Materials and Methods 

Plant material and DNA extraction 

Leaf material was obtained from different geograph- 
ical regions (Table 1). Figure 1 shows collection sites for 
those species whose preliminary phylogeographical pattern 



Table 1 - List of the species studied, their taxonomic position according to the infrageneric classification of Feuillet and MacDougal (2003), and 
GenBank accession numbers for DNA sequences. The number of individuals examined, the number of collection sites and the geographic origin of the 
populations are given for each species. 



Subgenus 


Species 


N 


Collection sites" 


Geographic regions'' 


Sequence types'^ 


Genbank 


Passiflora 


P. amettiystina 


3 


3 


RS, SC, SP 


3 


EU258307 - EU258309 




P. caerulea 


24 


14 


RS 


4 


EU258310-EU258313 




P. cinccinata 


7 


6 


MS, MG, BA, PE 


7 


EU258352 - EU258358 




P. edmundoi 


5 


4 


BA, MG, RJ 


5 


EU258370 - EU258374 




P. edulis 


12 


8 


RS, SC, MG 


10 


EU258375 - EU258384 




P. eiciileriana 


7 


4 


RS. SC 


3 


EU258317 -EU258319 




P. foetida 


9 


6 


RS, SC, PE, PB 


7 


EU258386-EU258391; 
EU258394 




P.jilekii 


5 


4 


SC 


1 


EU258320 




P. miersii 


3 


1 


SP 


1 


EU258322 




P. sidaefolia 


12 


2 


SP, MG 


1 


EU258435 




P. tenuifila 


15 


2 


RS 


8 


EU258446 - EU258453 




P. uruhiciensis 


12 


11 


RS, SC 


1 


EU258326 




P. villosa 


6 


5 


SC, PR, SP, MG 


6 


EU258392 - EU258393; 
EU258466 - EU258469 


Decaloba 


P. capsularis 


42 


5 


RS, SC, PR, SP, MG 


19 


EU258327 - EU258345 




P. misera 


15 


5 


RS 


4 


EU258409 - EU258412 




P. morifolia 


7 


5 


RS, SP 


2 


EU258323 - EU258324 




P. organensis 


17 


9 


RS, SC, PR, MG 


13 


EU258414-EU258426 




P. pohlii 


5 


4 


MG, MS 


1 


EU258325 




P. tricuspis 


8 


4 


SP, MS, DP 


6 


EU258455 - EU258460 




P. urnaefolia 


9 


4 


RS, SP, MG 


4 


EU258462 - EU258465 


Astrophea 


P. haematostigma 


19 


6 


SC, PR, SP, MG 


9 


EU258395 - EU258403 




P. rhamnifolia 


12 


5 


MG 


6 


EU258427 - EU258429; 
EU258432 - EU258434 


Deidamioides 


P. ovalis 


18 


7 


ES, BA 


11 


EU258359 - EU258369 



The collection sites for P. cinccinata P. organensis and P. haematostigma are shown in Figure 1 . 

'Abbreviations refer to Brazilian states: RS: Rio Grande do Sul; SC: Santa Catarina; PR: Parana; SP: Sao Paulo; RJ: Rio de Janeiro; ES: Espirito Santo; 
MG: Minas Gerais; MS: Mato Grosso do Sul; DF: Distrito Federal; BA: Bahia; PE: Pemambuco. 
"Sequence types found for the whole region (ITSl + 5.8S + 1TS2). 
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Figure 1 - Map of Brazil indicating the origin of samples of the tliree Passiflora species whose preUminaiy phylogeographical patterns were evaluated. 
Squares - P. cincinnata, circles - P. organensis, and triangles - P. haematostigma. 



was evaluated. Voucher specimens were deposited in the 
ICN Herbarium, Botany Department, Biosciences Institute, 
Federal University of Rio Grande do Sul. Total DNA was 
extracted from young leaves dried in silica gel, using the 
method of Roy et al. (1992). 

PGR amplification and sequencing 

Internal transcribed spacers (ITS 1 and 2) were ampli- 
fied using primers and amplification conditions as de- 
scribed by Desfeux and Lejeune (1996). To exclude the 
presence of low stability templates, 1 0% dimethyl sulfo- 
xide (DMSO) was used (Buckler IV et al, 1997; Fuer- 
tes-Aguilar et al, 1999). PGR products were checked by 
horizontal electrophoresis in 1% agarose gel stained with 
GelRed (Biotium). All PGR products were purified using 
the polyethyleneglycol (PEG) precipitation method (Dunn 
and Blattner, 1987). The sequencing reaction was per- 
formed as described by Sanger et al{\911) using ET termi- 
nators kit (GE Healthcare) on a MegaBAGE 1000 auto- 
matic sequencer (Amersham Biosciences). The quality of 
the sequences was examined through the Ghromas package 
(available from http://www.technelysium.com.au/ 
chromas.html). Nucleotide sequences were searched 
against the GenBank database (http://ncbi.nlm.nih.gov/ 
BLAST) using BLAST tools (Altschul et al, 1990). 

Data analysis 

GenBank numbers for the obtained sequences are 
given in Table 1 . Forward and reverse reads were analyzed 
for all sequences. A site was identified as ambiguous when 



double peaks occurred in the same position in both strands, 
with the weakest signal reaching at least 25% of the 
strength of the strongest (Fuertes-Aguilar et al, 1999; 
Fuertes-Aguilar and Nieto-Feliner, 2003). The term 'am- 
biguous' was used instead of 'heterozygote site', as the ori- 
gin of ambiguity could not be determined with any degree 
of certainty, since a variation might represent an actual 
heterozygote (for which each allele had been inherited from 
one parent) or only a variation among different ITS copies 
across the genome. Thus, in all analysis, ambiguous sites 
were treated as missing data. Sequences were aligned with 
the PRANK program (Loytynoja and Goldman, 2008). 
Variable sites, nucleotide and haplotypic diversity, were 
estimated using the Arlequin 3.11 (Excoffier et al, 2005) 
software. Relationships among sequences were inferred by 
means of median-joining networks (Bandelt et al, 1999) 
under Network version 4.5. Finally, inter and intraspecific 
distances were estimated with PAUP* version 4.0 
(Swofford, 1998), by taking the distance estimated through 
maximum likelihood, and assuming the evolutionary 
model suggested under the Aikaike Information Criterion 
in the Modeltest 3.7 program (Posada and Grandall, 1998). 
Since the sequences among subgenera were too divergent 
to be combined in the same alignment, diversity in each 
subgenus was analyzed separately. 

Results 

GC content ranged from 49% to 66% (Table 2 - align- 
ments are available directly from the authors), with an aver- 
age of 60%, being significantly higher in the subgenus 
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Passiflora (63%) than in Decaloba (53%) (p < 0.002). ITS 1 
length varied between 225 bp {P. eichleriana Mast.) and 
278 bp {P. misera HBK), whereas for ITS2, this variation 
was from 1 84 bp {P. sidae folia M. Roemer) to 222 bp {P. 
caerulea L.) (Table 2). ITS length was significantly differ- 
ent between the two subgenera for both ITS 1 and ITS2 (Ta- 
ble 3). The size of the 5.8S gene (160 bp) was constant in all 
analyzed species, although a certain intraspecific polymor- 
phism, half of which represented by ambiguous sites, was 
found in this region in P. capsularis L., P. edulis Sims., P. 
urnaefolia Rusby. and P. villosa Veil. 

With the exception of P. pohlii Mast., in Mart, (sub- 
genus Decaloba) and P. urubiciensis Cervi (subgenus 
Passiflora), both of which were monomorphic, all species 
had some degree of intraspecific variability, represented by 
insertion/deletion (indels) and substitution events, part of 
which involving ambiguities. When we excluded all ambi- 
guities from the dataset for subsequent analyses we found 
that three more species become monomorphic {P. jilekii 
Wawra, P. miersii Mast, in Mart, and P. sidaefolia, all be- 
longing to the subgenus Passiflora). The number of nucleo- 
tide substitutions was higher than indels, and transitions 
more numerous than transversions (Table 2). The intra- 
specific index of nucleotide diversity varied from zero to 
0.0227 (Table 2). The highest intraspecific pairwise dis- 
tance was 0.0447 {P.foetida L.). Pairwise interspecific dis- 
tances between individuals belonging to species from the 
same subgenus, varied from 0.0128 to 0.0260 in the subge- 
nus Astrophea (DC.) Mast, (both between P. 
haemato stigma Mart, ex Mast, and P. rhamnifolia Mast.), 
from 0.0049 {P. tricuspis Mast, in Mart. vs. P. organensis 
Gard.) to 0.3860 {P. capsularis vs. P. morifijlia Mast, in 
Mart.) in Decaloba and from zero {P. miersii Mast, in Mart. 
vs. P. edmundoi Sacco) to 0.1713 {P. caerulea vs. P. 
ftyetida) in Passiflora. For the subgenus Deidamioides 
(Harms) Killip a single species was studied. The distribu- 
tions of inter and intraspecific pairwise distances are shown 
in Figure 2. 

As in this study we presented an overview of ITS ge- 
netic diversity in Passiflora species as a whole, rather than 
an exhaustive analysis of just a few species, it could be ar- 
gued that the somewhat limited, and to some extent, un- 
equal sampling of the different species may tend to bias our 
estimates. To assess the influence of sample size on genetic 
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Figure 2 - Histogram showing the distribution of frequency values of intra 
(empty bars) and interspecific (solid bars) genetic distances. Interspecific 
genetic distances were calculated only comparing species belonging to the 
same subgenus. 



diversity, we calculated the Spearman correlation coeffi- 
cient between either sample size or number of sampled 
populations against two measures of genetic diversity, the 
number of different sequence types and nucleotide diver- 
sity. We observed a significant positive correlation be- 
tween sample size and both the number of sampled popula- 
tions and that of different sequence types (rS = 0.601, 
p = 0.002 and rS = 0.462, p = 0.027, respectively). Never- 
theless, when nucleotide diversity was considered, no sig- 
nificant correlation was found with either sample size or the 
number of sampled localities (rS = 0.070, p = 0.751; and 
rS = 0.260, p = 0.231; respectively). These results indicate 
that the values presented here for nucleotide diversity are 
independent of sampling artifacts, and thus may be viewed 
as approximations of the level of genetic diversity in these 
species. However, it is also evident from our data that in- 
creased sampling would be required for a better character- 
ization of all sequence types present in a given species. 

When considered together with sampling locations, 
genetic variation for at least six species {P. 
haematostigma, P. organensis, P. cincinnata Mast., P. 
edmundoi, P. villosa, and P. capsularis) showed some 
geographic structure as can be observed from the network 
connecting the different sequence types (Table 4; Figu- 
re 3). For P. cincinnata (Figure 3a) there was a central 
core of sequence types located in NE Brazil (BA and PE 
states), to which the others were connected, the most dis- 
tantly located population also having the most divergent 
sequence type. A more complex picture was observed for 
P. organensis (Figure 3b), wherein the most distinct lin- 



Table 3 - Passiflora subgenera average characteristics. Values within parentheses refer to the range of variation or standard deviation. 



Subgenus 


Species number 


%GC 


ITSl length 


ITS2 length 


Ambiguous sites 


Nucleotide 


Sequence 












number 


diversity 


diversity 


Passiflora 


13 


63* (60-66) 


227* (225-233) 


216* (184-222) 


4 (0-22) 


0.0062 (0.0036) 


0.5489 (0.0611) 


Decaloba 


7 


53* (49-54) 


276* (271-278) 


200* (198-203) 


3 (0-8) 


0.0021 (0.0014) 


0.5736 (0.0595) 


Astrophea 


2 


63 


269 


206 (204-208) 


5 (3-7) 


0.0035 (0.0022) 


0.7842 (0.0564) 


Deidamioides 


1 


64 


236 


212 


7 


0.0048 (0.0028) 


0.8635 (0.0433) 



* Statistically significant values (p < 0.002). 



104 



Mader et al. 



Table 4 - Collection points of Passiflora species considered in phylogeographic analyses. 



Species Sampling sites Geographic coordinates Sequence types 



P. cincinnata 


1- CasaNova/BA 


09° 09'43" S/40° 58'15" W 


c2 




2- Petrolina/PE 


09° 23'55" S/40° 30'03" W 


cl 




3- Abaira/BA 


13° 14'59" S/41° 39'49"W 


c3,c4 




4- Caetite/BA 


14° 07'36" S/42° 26'05" W 


c3 




5- Brumadinho/MG 


20° 08'36" S/44° 11'59"W 


c7 




6- Jardim/MS 


21° 28'49" S/56° 07'53" W 


c6 


P. organensis 


7- Caete/MG 


19° 52'56" S/43°40'12"W 


o3 




8- Munhoz/MG 


22° 36'31" S/46° 2r4r' W 


ol, o2 




9- Quatro Barras/PR 


25° 20' 16" S/48° 54'49" W 


o4, o5, o6 




10- Petrolandia/SC 


27° 35'29" S/49° 44' 14" W 


o5, o8, o9 




1 1 -Bom Jardim da Serra/SC 


28° 22'14" S/49° 34'41" W 


o8 




12- Timbe do Sul/SC 


28° 49'42" S/49° 50'30" W 


o8 




13- Praia Grande/SC 


29° 07'25" S/49° 58'29" W 


o7 




14- Tres Cachoeiras/RS 


29° 26'53" S/49° 55'37" W 


olO 




15- Maquine/RS 


29° 34'36" S/50° 17'35"W 


oll,ol2,ol3 


P. haematostigina 


16- Sao Gonpalo do Rio Preto/MG 


18° 00' 15" S/43°23'27"W 


h2,li4 




17- Araiijos/MG 


19° 53'59" S/45° 15'2r' W 


h2 




18- Itabirito/MG 


20° 22'01" S/43° 39'40" W 


hl,h2, h3,h4 




19- Biritiba/SP 


23°39'04" S/46° 08'17"W 


h5, h6, h7 




20- Guaratuba/PR 


25° 52'58" S/48° 34'29"W 


h9 




2 1 - Rancho Queimado/SC 


27° 40'59" S/49°01'59"W 


Ii8,li9 




Figure 3 - Median-joining networks based on ITS sequence types. Each circle represents a sequence type, their sizes being proportional to respective fre- 
quencies. Shading indicates the state in Brazil where each sequence-type was found (according to boxes; key to the abbreviations in Table 4). Transversal 
bars indicate the number of mutations that differentiate sequence types; Mv: median vector. The differences between the numbers of sequence types indi- 
cated in Table 1 for P. haematostigina are due to hypervariable sites that were removed for network analysis. 3a -P. cincinnata; 3b - P. organensis; 3c -P. 
haematostigma. The relationship between sequence types and collection sites is indicated in Table 4. 



eages were found at the two extremes of distribution (MG 
and RS, see Figure 1), whereas another group of lineages 
was arranged in a roughly north/south gradient. A clearer 
pattern emerged for P. haematostigma (Figure 3c), with a 
clear separation between the sequences from Minas Ge- 



rais (MG) and those from other more southern localities, 
which are arranged in a north/south gradient (SP, PR, SC). 
Finally, in the three remaining species we also found dif- 
ferences among isolated individuals from different re- 
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gions albeit with a non-trivial relationship between 
genetic lineages and geography. 

Discussion 

In almost all the species sampled in this study (21 out 
of 23), a certain degree of genetic variation occurred in the 
ITS region, even though for some of these only few individ- 
uals were sequenced. This is in agreement with other stud- 
ies, thereby implying that in Passiflora, the ITS region 
seems to be more informative at the populational level than 
other marker. Lorenz-Lemke et al. (2005), when studying 
32 individuals of P. actinia Hook and 20 individuals of P. 
elegans Mast., and Koehler-Santos et al. (2006), studying 
32 plants from P. alata Curtis, all discovered genetic varia- 
tion on using ITS but none with cpDNA markers. Lorenz- 
Lemke et al. (2005) found 32 polymorphic sites. Nucleo- 
tide diversity was 0.0060 in P. actinia and 0.0020 in P. 
elegans, whereas Koehler-Santos et al. (2006) detected an 
ITS nucleotide diversity value of 0.0036 as in P. alata. 

Patterns of intraspecific variation in Passiflora have 
been also studied by using different markers. More specifi- 
cally, Fajardo et al. (1998) employed RAPD markers and 
found higher levels of variation in P. ligularis Juss. and P. 
adenopoda DC. when compared to P. edulis and P. 
maliformis L., whereas Sanchez et al. (1999) used cpDNA 
restriction fragment length polymorphism, and found intra- 
specific variation in four species (P. maliformis, P. 
ligularis, P. edulis, and P. mollissima (Kunth) Bailey), but 
no variation in another three {P. caerulea, P. sp. India, and 
P. adenopoda). Finally, a morphometric approach under- 
taken by Plotze et al. (2005) found different levels of varia- 
tion in leaf vein pattern among ten species, with P. caerulea 
being the most diverse. Even though differences among 
sampled species and markers make a direct comparison of 
these studies difficult, the emerging overall picture is that 
intraspecific variation is not evenly distributed among spe- 
cies, and that a common pattern of intraspecific diversity 
for the nuclear genome, the plastid genome and morphol- 
ogy may be hard to obtain, given the discrepancy of differ- 
ent datasets for the same species {e.g. P. caerulea, P. edulis 
and P. maliformis'). Possibly, this picture reflects the com- 
plexities of the evolutionary history of the genus, and indi- 
cates that robust patterns would only emerge when 
different markers are considered together. 

In this study we did not take into consideration sites 
for which two nucleotides could be detected in the same in- 
dividual. Because there are several ITS repeats and because 
each individual inherited its genome from both of its par- 
ents, it is not possible to know whether sites harboring two 
different nucleotides are really heterozygous or whether 
they rather represent ITS sequence variation across differ- 
ent repeats in a single haploid genome. It is therefore im- 
possible to infer meaningful haplotypes based on ITS 
sequences obtained from PCR amplifications from total 
DNA. Importantly, because we excluded from our analysis 



all sites with ambiguous information, the degree of genetic 
variation in the ITS region in this genus is actually an un- 
derestimate. A way of maximizing the detection of genetic 
variation in ITS sequences would be by cloning all individ- 
uals and sequencing multiple clones for each. However, 
this alternative is not only expensive for phylogeographic 
and population genetics studies, but might also be subject to 
criticism because one would be mixing orthologous and 
paralogous variation in a single analysis (Bailey et al., 

2003) . AK Hansen (PhD. Thesis, University of Texas, 

2004) sequenced ITS clones obtained for 41 Passiflora spe- 
cies and observed that in 25% of the studied taxa, various 
ITS sequences were more related to sequence types in other 
species. It is not known whether this illustrates the com- 
plexity of using paralogue ITS sequences in a single analy- 
sis or not, since, according to other studies on the same 
genus, variation patterns inconsistent with sequence types 
from a single monophyletic species have never been found 
(Muschner et al, 2003; Krosnick and Freudenstein, 2005; 
Lorenz-Lemke et al, 2005). 

A comparison of intra and interspecific genetic dis- 
tances (Figure 2) suggests little overlap between the two 
levels, with more than 50% of the intraspecific compari- 
sons resulting in values larger than zero, thus stressing the 
potential of the ITS region as a molecular marker in phylo- 
geographic studies in Passiflora. The narrow overlap be- 
tween intra and interspecific values may be understood in 
two-ways. First, the intraspecific genetic variation may be 
inflated if there are unrecognized species within a given 
taxon name. This may be the case for P. foetida, for which 
38 varieties are recognized (Killip, 1938), and which 
showed the highest intraspecific variation (Table 2). A 
thorough taxonomic revision of this species would there- 
fore considerably reduce the overlap between intra and 
interspecific genetic variation in the ITS region. Second, 
interspecific distances are expected to be low between re- 
cently diverged species. Thus, we may expect that either re- 
cent speciation events or taxonomic confusion would cause 
an overlap between genetic distances estimated at these two 
taxonomic levels. This may be of special concern for 
speciose groups. Moreover, differences represented by 
indel events are not taken into account when estimating dis- 
tances, thus probably reducing the overall estimate of inter- 
specific divergence. For example, some individuals fromP. 
miersii and P. edmundoi are only differentiated by a 2-bp 
indel in the ITS2 region. 

Provided that concerted evolution is often invoked as 
an explanation for the maintenance of genetic homogeneity 
among ITS copies, why then does genetic variation in the 
ITS region seems to be the rule rather than the exception in 
Passiflora species? The time for homogenizing different 
ITS copies depends on the number of the different copies, 
on the number of different chromosomes containing 
nrDNA, and on the crossing-over rate among chromosomes 
(Quijada et al, 1998), since crossing-over towards telo- 
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meres seems to be more common, as previously shown in 
Gossypium L. and Thinopyrum A. Love (Wendel et al., 
1995; Li and Zhang, 2002). Melo and Guerra (2003) 
mapped the 45 S ribosomal DNA (including the ITS region) 
in several Passiflora species. In none of the species ana- 
lyzed (including P. capsularis, P. cincinnata, P. edulis, P. 
haematostigma and P. morifolia,) had nrDNA located in 
telomeres, this indicating that chromosome location in this 
genus may be, at least partly, responsible for a slow down in 
the homogenization of the different ITS copies. This could 
partially explain, independent of the taxonomy issues, the 
high values observed for P. foetida, which has six nrDNA 
clusters far from telomeres, in regions of reduced cross- 
ing-over. 

The differences in ITSl and ITS2 length and GC con- 
tent, detected among species in the subgenera Passiflora 
and Decaloba, might be related to their significantly differ- 
ent evolutionary rates, as shown by Muschner et al. (2003). 
The latter also ascertained that ITS 1 and ITS2 nucleotide 
diversities were significantly higher in Decaloba than in 
Passiflora species. Decaloba species are in general self- 
compatible (Ulmer and MacDougal, 2004), a condition that 
favors lower intrapopulation and higher interpopulation ge- 
netic diversities, as there is less pollen dispersion (Maki et 
al, 1999). In the present work, no difference was detected 
when average intraspecific nucleotidic diversities were 
compared in Passiflora. This indicates that, despite the dif- 
ferent evolutive rates, ITS intraspecific variability it is not 
distinct in these subgenera. We consider that the differenti- 
ated evolutionary rates in the subgenera Decaloba and 
Passiflora do not represent obstacles to the use of ITS 
markers for intraspecific populational analysis in 
Passiflora species. However, this question must be handled 
carefully when dealing with phylogenetic questions in this 
genus and is worthy of future study. 

Even though our study did not aim at detecting fine 
phylogeographic structure for these species, for at least six 
species, sampling site location could be related to the ob- 
served genetic diversity. Despite the low number of sam- 
ples, the overall pattern seems to indicate a general 
latitudinal trend. Of course, detailed phylogeographic stud- 
ies will be required to reveal accurate geographic patterns, 
and to infer those evolutionary processes that influenced 
the distribution of the genetic lineages currently found in 
these species. It is interesting to note that Lorenz-Lemke et 
al. (2005), also based on ITS sequences, detected similar 
patterns in P. actinia and P. elegans. Palynological data in- 
dicate two Holocene expansions in the Brazilian Atlantic 
Rain Forest from the northeast towards south Brazil 
(Lorscheitter, 1997). Nonetheless, any extensive assess- 
ment of whether the history of the genetic lineages within 
these species indeed reflects concordant colonization 
routes in the Brazilian Atlantic Rain Forest will await fur- 
ther studies and additional data. 



Despite their complex evolution, ITS sequences have 
high inter and intraspecific variability in Passiflora species. 
This information can be of importance for accessing their 
evolutionary history and contributing to conservation. At 
least in some species, a clear geographic stmcture of se- 
quence types can be identified, with expressive agreement 
among the different species. As a whole, the results pre- 
sented here indicate that ITS is a useful tool for evaluating 
intraspecific variability in Passiflora species. 
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